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A quasi-potential approximation to the Navier-Stokes equation for low viscosity fluids 
is developed to study pattern formation in parametric surface waves driven by a force 
that has two frequency components. A bicritical line separating regions of instability 
to either one of the driving frequencies is explicitly obtained, and compared with ex¬ 
periments involving a frequency ratio of 1/2. The procedure for deriving standing wave 
amplitude equations valid near onset is outlined for an arbitrary frequency ratio follow¬ 
ing a multiscale asymptotic expansion of the quasi-potential equations. Explicit results 
are presented for subharmonic response to a driving force of frequency ratio 1/2, and 
used to study pattern selection. Even though quadratic terms are prohibited in this case, 
hexagonal or triangular patterns are found to be stable in a relatively large parameter 
region, a fact that is in qualitative agreement with experimental results. 


1. Introduction 

When a fluid layer is periodically oscillated in the direction normal to the free surface 
at rest, parametric surface waves (or Faraday waves) appear above a certain critical 
value of the vibration amplitude (Miles & Henderson 1990; Cross & Hohenberg 1993). 
We present in this paper an extension of a weakly nonlinear model previously introduced 
by Zhang & Vinals (1996a) and (1996b), which is valid in the limit of low fluid viscosity, 
large aspect ratio and large depth, to study pattern selection near onset of Faraday waves 
when the driving force has two independent frequency components. 

The model developed by Zhang & Vinals (1996b) was based on a quasi-potential ap¬ 
proximation to the equations governing fluid motion. In it, the flow is considered to 
be potential in the bulk, subject to effective boundary conditions at the moving free 
surface that incorporate the effect of the rotational component of the flow within a thin 
boundary layer near the free surface. We further assumed without rigorous justification 
that for low viscosity fluids, only linear viscous terms need to be retained in the resulting 
equations (the so-called “linear damping quasi-potential equations”, or LDQPE’s). A 
multi-scale analysis of the resulting LDQPE’s led to the prediction of standing wave pat¬ 
terns of square symmetry near onset for capillary waves, in agreement with experiments. 
For mixed capillary-gravity waves, patterns of hexagonal symmetry or quasi-periodic 
patterns were predicted depending on the value of the damping coefficient. 
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Although some of our predictions have been confirmed experimentally (Kudrolli & 
Gollub 1996), we address in this paper pattern selection in systems driven by periodic 
forces comprising two frequency components for two reasons. First, recent detailed ex¬ 
periments involving two frequencies provide a good opportunity for additional tests of 
our theory. Second, additional control parameters appear relative to the single frequency 
case, namely, the frequency ratio, relative amplitude and phase difference. Therefore, 
one may anticipate richer dynamics and more interesting steady states as compared to 
the single frequency case. 

Consider a fluid layer perpendicular to the 2 : axis, driven by a force of the form 


gz{t) = -go - gz[r sm2muJot + (1 - r)sin{2nuJot + (j))], (1.1) 


where r (0 < r < 1) is the relative amplitude of the two frequency components, 4> is their 
relative phase difference, go is the constant background gravitational field pointing to the 
negative 2 : direction, and is the amplitude of the oscillatory component. When r is 
sufficiently small, the forcing component sin(2na;ot-l-((') is expected to dominate, the other 
component (sin2mwot) being a small perturbation. Linear stability of the flat surface 
is thus approximately determined by subharmonic instability to the forcing component 
sin(2na;ot + ^) (note that this is not necessarily a subharmonic response to the entire 
driving force), with a critical wavenumber kon given by gokon + F/cq^/p = where 

F is the surface tension of the surface and p is the density of the fluid. On the other 
hand, when r is sufficiently close to 1, the forcing component sin2mtuot dominates, and 
linear stability is approximately determined by subharmonic instability at a frequency 
muo (again it is not necessary that it be a subharmonic response to the entire driving 
force), with a critical wavenumber kom given by gokom + = tu^Wq. When r is 

varied from 0 to 1, one could imagine that the above two instabilities can have the same 
threshold value of the driving amplitude gz for some value of r, which will be referred to 
as the hicritical value rt. Furthermore, given a fixed frequency ratio mjn, the bicritical 
value Tb can be a function of the phase difference </>. In the r — (j) plane, the bicritical 
values rb{(j)) form a line (the bicritical line) which separates the two regions with different 
characteristic temporal dependencies and spatial scales. 

The first two-frequency Faraday experiments were reported by Edwards and Fauve 
(Edwards & Fauve 1992, 1993 and 1994), albeit in fluids of large viscosity. In that 
case, and with a purely sinusoidal parametric force, roll patterns (or lines) similar to 
those observed in Rayleigh-Benard convection are observed near the primary instability 
of a flat surface. In the case of two frequencies, Edwards and Fauve studied several 
ratios m/n, including 3/5, 4/5, 4/7, 6/7, and 8/9, with most of their data for the 
case m/n = 4/5. In this latter case, parametric surface waves were found to respond 
synchronously (harmonically) with the force when the even frequency forcing component 
sin(2 X dwot) dominates, and subharmonically when the odd frequency forcing component 
sin(2 X 5uJot + 4') dominates. The bicritical line obtained by Edwards and Fauve is almost 
independent of the phase difference </>, and for m/n = 4/5, rb{(j)) ~ 0.32. 

They also found many interesting standing wave patterns near the primary instability 
of a planar surface: lines, squares, hexagons, and twelve-fold quasi-crystalline patterns. 
For the cases of m/n = even/odd (4/5, 4/7, 6/7, and 8/9), hexagonal patterns were 
observed for a wide range of values of r, in the region of harmonic response in the r — (/> 
plane (when the even frequency forcing component dominates). On the same side of the 
bicritical line where hexagonal patterns are found, a stable twelve-fold quasi-crystalline 
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pattern was observed in a small region very close to the bicritical line. For min = 4/5, 
this small region is around r 0.32 and </> w 7.5°. |t| 

Two-frequency driven Faraday experiments using less viscous fluids have been reported 
more recently by Muller (1993). In the case of single frequency forcing, Muller observed 
standing wave patterns of square symmetry near onset, in agreement with previous ex¬ 
periments (Lang 1992; Ezerskii, Rabinovich, Reutov & Starobinets 1986; Tufillaro & 
Gollub 1989; Ciliberto, Douady & Fauve 1991; Christiansen, Alstrp m & Levinsen 1992; 
Bosch & van de Water 1993; Edwards & Fauve 1993) and our previous theoretical work 
(Zhang & Vinals 1996b) for weak viscous dissipation. He also studied pattern formation 
near onset in a system driven by a two-frequency force, with frequency ratio m/n = 1 / 2 , 
and found that the r — ip plane is divided into two regions: in the region of larger r, para¬ 
metric surface waves are found to respond subharmonically to the entire driving force, 
whereas in the region of smaller r surface waves respond harmonically. The bicritical 
line has an interesting dependence on the phase difference which was not observed 
by Edwards and Fauve (1993) and (1994) for frequency ratios other than 1/2 in highly 
viscous fluids. 

In addition to square patterns, Muller found hexagonal patterns that respond harmon¬ 
ically to the driving force, and hexagonal or triangular patterns that respond subhar¬ 
monically to the driving force. Interestingly, a spatially disordered region, that responds 
subharmonically to the driving force with an additional slow time dependence, was also 
found near onset. 

We study in this paper how much of this seemingly complicated stability diagram found 
by Muller can be understood, at least qualitatively, in terms of perturbative analysis 
of the quasi-potential equations introduced by Zhang & Vinals (1996a) and (1996b). 
In Muller’s experiments, the damping parameter (7 = 2iyk^, where u is the kinematic 
viscosity of the fluid, and k the wavenumber of a surface Fourier mode) 7 s = 0.17 in the 
region of subharmonic response, and jh = 0.33 in the region of harmonic response. Since 
the quasi-potential approximation is only valid in the limit of small viscous dissipation 
(small 7 ), we anticipate qualitative agreement with experiments, while smaller damping 
parameters are needed to quantitatively test our theoretical predictions. 

The rest of this paper is organized as follows. In section we present analytical results 
for the bicritical lines, and make comparison with Muller’s experimental results. Section 
^ contains the derivation of standing wave amplitude equations (SWAE’s), and results of 
pattern selection based on the SWAE’s. Further discussion and conclusion are presented 
in the section 


2. Bicritical line 

In this section, we consider the linear stability of two-frequency forced Faraday waves, 
and obtain analytical results for the bicritical line rb{4>). We restrict our attention to 
cases in which m and n are two small positive integers that are relatively prime (e.g, 
m/n = 1/2,1/3, 2/3). Without loss of generality, we assume m < n. The phase difference 
(j) can be chosen within —tt/to < </> < tt/to because g(t) is invariant with respect to the 
transformation: ^ -1- 2 tt/to and t ^ t + jir/mujo for any integer j that makes 


f In the form of the driving 

force used by Edwards and Fauve, —go + ao[cos x cos(4t<;t) -|- sin y cos(5:^t -I- the twelve-fold 
quasi-crystal line pattern w as observed around y « 65° and (j)' ~ 75°. Also, for later reference, 
we note that Muller (199311 used the definit ion a{t\ = — go-l-flz[r cos 2aJot-l-(1 — r) cos(4aTot-l-(p)]; 
which can be written in tne form of Eq. (1.1) by the transformation oi t ^ t + 7r/(4a;o) and 
p = 4> — 7r/2. 
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{nj + l)/m an integer. Linear stability of the free surface of a weakly viscous fluid 
of infinite depth under a two-frequency parametric force is determined by the damped 
Mathieu equation (Landau & Lifshitz (1976)), 


dtthk + ^vk^dthk + 


Tk^ 

gok-\ - \-gzk{ rsm2mLUot+(l — r) sm{2nuJot + (j)) 

P 


hk=Q. ( 2 . 1 ) 


2.1. Subharmonic versus synchronous responses 

The two-frequency parametric force in Eq. (H has an angular frequency 2loo, or a 
period tt/uiq (recall that m and n are assumed to be relatively prime). Whether a 
region is subharmonic or synchronous (harmonic) with respect to the total driving force 
depends on whether m/n is odd/odd, even/odd, or odd/even. This overall response to 
the total force is important since it can be easily determined experimentally by using 
stroboscopic methods. According to Floquet’s theorem, the general form of the solutions 
of Eq. ( |2.1| ) is a periodic function multiplied by an exponential function of time. At 
the stability boundaries, Eq. (2.1) has periodic solutions, hk{t + tt/wq) = phk(t), where 
the Floquet multiplier p = ±1. For p = —1, hk{t) has a period of 27r/a;o, and thus 
the response is subharmonic, whereas for p = 1, hk{t) has a period of tt/wo, and hence 
the response is harmonic. Therefore the solution of Eq. (I3) at the boundaries for 
subharmonic instability {p = —1) can be written as a Fourier series in odd frequencies 
(2j — l)ujQ,j = 1,2 ,..., oo, whereas at the boundaries for harmonic instability (p = 1) 
the frequencies involved in the series are 2ja;o, j = 1, 2 ,..., oo. As a consequence, when 
the odd (even) frequency dominates, the response is subharmonic (harmonic). 

We now turn to a detailed calculation of the bicritical line in the (r, (p) plane, starting 
from the linearized quasi-potential equation (Eq. (p^)) (Zhang & Vinals 1996a and 
1996b). Since the characteristic time and length scales are different on the two sides 
of the bicritical line for two-frequency driven Faraday waves, we shall use dimensional 
variables, and group them into dimensionless quantities when necessary. In ord er to keep 
our notation simple, we define 5 = g^k J-Tkr’/p, and / = gzk/A. Equation ( 0 ) can now 
be written as 


dtthk + ‘^jdthk -l- 


6 + Af I r sin 2 mi 0 ot -|- (1 — r) sin(2na;ot -I- </>) 


hk = 0 . 


( 2 . 2 ) 


2.2. Multiple scale expansion 

We consider small values of the driving amplitude / and of the damping coefficient 7 , 
and introduce a small parameter p such that 7 = p 7 o and / = p/o- For simplicity, 
we will consider the wavenumber /cq exactly at subharmonic resonance to either of the 
forcing components, and thus no expansion for the wavenumber is needed, define 60 = 
goko -|- Efcg/p, and seek a solution hk in a power series in p as, 

hk = h<f\t, Ti, T 2 ) + phi^\t, Ti, T 2 ) + v^h<'^\t, Ti, T 2 ) + • • •, (2.3) 

where Ti = rjt and T 2 = rf't. We have introduced a second slow time scale T 2 in the 
expansion for hk, in addition to the slow time scale Ti. The second slow time scale is 
necessary since we are going to perform our perturbation expansion up to 0 {rf'), which 
is one order higher in p than the perturbation expansion for the case of a single sinusoidal 
driving force in ^hang fc Vinals (1996b) . 

Equation ( | 2 . 2 D contains a damping term proportional to 7 , while terms of 0 ( 7 ^/^) and 
0 ( 7 ^) (or 0(jf)) have been neglected in the quasi-potential approximation (Zhang & 
Vinals 1996b). It seems to be inconsistent to consider terms of 0(p^) in the solutions to 
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Eq. ( |2.2| ). The basic reason to consider solutions up to 0{rf) here is that, as we shall 
see below, these terms can have relatively large coefficients due to the special nature of 
the expansion. When the dominant driving frequency is 2 towo, the small parameter in 
this perturbation expansion is the dimensionless driving amplitude r//(m^WQ) 1. The 

driving force in this case can be written as 


4r/ 

9 2 


1 — r 

sin2mwot H- sin{2nuJot + (p) 


(2.4) 


Since the perturbative expansion is in be in powers of rf / (to^Wq), as well as in (1 — r)/r, 
it is necessary to assume that (1 — r)/r ~ 0{1) or smaller. However, (1 — r)/r can be quite 
large on the bicritical line if r = ri,{(j)) is small for some values of p. In this case, higher 
order terms in rj can be important if they have factors of (1 — r)/r in their coefficients. 
In Muller’s experiments, for example, the smallest value of r{,((/)) is about 0.2, and thus 
(1 — r)/r ~ 4. Similar arguments can be made when 2nu)Q is the dominating driving 
frequency. In the case of Muller’s experiments, the largest value of r/(l — r) along the 
bicritical line is less than 0.5. A consistent calculation would require keeping terms of 
0 ( 7 ^/^) or higher at the linear level of approximation in the surface variables (Eq. dH), 
leading to a much more involved nonlinear analysis. However, the agreement that we 
find with Muller’s bicritical line is quite reasonable. 

On substituting the expansion for hk into Eq. (2.2), we have at 0{rp), 


dtthf + Sohf = 0 . 


(2.5) 


At 0 ( 77 ), we have, 




HO) 


“ dh 

-4/o( rsin2TOWot+ (1 - r) sin{2nu;ot + (p) ] 


dTidt 


dt 


and at 0 {t]'^), we have. 


k ’ 

( 2 . 6 ) 




1(2) , A 1,(2) '2.^. u n^.' ^ dh 


( 0 ) 


dT2dt dTf 


- 270 - 


-2 




dTi dTidt 


- 270 - 


dtth):> + 6oh):> = -2 


We now consider separately the two cases of So = m^cog and Jq = n^ojo- 


dh 


( 1 ) 


dt 


— 4/0 ( r sin 2 niuJot + (1 — r) sin(2naJof + (p) ) . 


(2.7) 


3o = 

At 0 ( 77 °), we have the solution h® = A(Ti, T 2 ) cosTTiwot + i3(Ti, T 2 ) sinTuwot. At 
0 ( 77 ), we have the following solvability condition in order to avoid secular terms, 

dA ^ ^ r/o 
dTi ~ 




rrujjo 


dB _ _ / r/o 


dTi 


The solution at this order reads. 


- - 
"-k — 


\mijJo 

rfo 


+ I 0 ] B. 


( 2 . 8 ) 

(2.9) 


477i2a;Q 


(B cos 3muJot — A sin imujot) 


(1 - r)fo 
2 n{n — m)ujg 


B cos ( {2n — m)uJot + p] + A sin ( {2n — m)ujot + p 
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(1 - r)fo 
2n{m + n)ujQ 


B cos ( (2n + m)uJot + </>) — ^ sin ( (2n + m)uiQt + (j) 


( 2 . 10 ) 


Whether higher order contributions are important or negligible is not known at this 
point. Since we have found nontrivial equations for the amplitude A and B, and these 
equations give us a threshold value of the driving amplitude /o = mujo^ofr for the 
Faraday instability, one would guess that it is not necessary to consider higher order 
contributions. However, as we show below, contributions at 0{rf‘) are important to 
determine the bicritical line rb{4>). Interestingly, and for a quite obvious reason to be 
also discussed below, contributions at 0{rf) for the case of m/n = 1/2 are qualitatively 
different from the cases of any other frequency ratios. 

At 0{rf‘), we have. 


dtth^k^ + = 


2muJo 
2 m Wo 


dA 

m~ 

dB 

m 


d^B dB 


d^A 


■ 270 


dA 

HTi 


fU 2{l-rf 

Wq \ 2m? r? — 

fU 2{l-rf 


B 


sin mujfjt 


\2m^ 


A 


cos mujot 


2m? — mn + v? r(l — r)// 


2 nm^(n — 

m) 

V 

2 m^ + mn - 

F n^ r 

(l-r)/oY 

2 nm^(n + 

m) 

^ 0 ^ V 

2 (n^ + mn 

— m^) 

^’(1 - O/o 

mn{r? — 

m?) 

wg 

2{r? — mn 

— m?) 

^’(1 - r)f^ 


mn{r? — m^) 


° ' Acos[(2n — 3m)wot+(/] — Bsm[{2n — 3m)uJot+(j)] 
Acos[(2n + 3m)wot+(/)] + Bsm[{2n + 3m)uJot+(j)] 
A cos[(2n — m)ujot+(j)] + B sin[(2n — m)wot+(/>] 
A cos[(2n + m)wot+(/>] + B sin[(2n + m)wot+(/>] 


n(n — m)ujQ 

(l-r)Vg 

n(m + n)wg 


A cos[(4n — m)uJot + 2(j)] — B sin[(4n — m)uJot + 2(j)] 


A cos[(4n + m)ujot + 2(j>] + B sin[(4n + m)ujQt + 2(j>] 


-5—5- ( Afcos Smujot — cos 5 muJot] — Hfsin Smujot + sin Smwot] 

2m‘‘ujn V 


( 2 . 11 ) 


Terms on the RHS of the above equation will contribute to the solvability condition only 
if their frequencies are mwp. Thus, besides the first two terms, only the third term on 
the RHS of Eq. ( 2. Ill) contributes to the solvability condition when 2n — 3m = m, i.e., 
the frequency ratio m/n = 1/2 (note n > m by assumption). Therefore the solvability 
conditions at 0{rf') for frequency ratio m/n = 1/2 is different from that for all other 
frequency ratios. 

For m/n = 1/2, the solvability conditions read. 


dA 

m 

dB 

m ~~ 


1 

2wo 

1 

2wo 


r r 2 


3r^ 

3r^ 


2(1 — r)^ 


+ r(l — r) cos (/ ) — 7o 




2Wg 


( 2 . 12 ) 


2(1 — r)^ 


+ r(l — r) cos (/ ) — 7o 




2 wq 


(2.13) 
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For a frequency ratio mfn 1/2, the solvability conditions are, 


7 


dA 

1 

'III 

^ 3r^ 

W2 ~ 

2muJo 

.^0 ' 


dB 

1 

'1 

/ 3r^ 

W2 ~ ' 

2muJo 


\2m^ 


2(1 

2 (l-r)^ \ 

— w? ) 



(2.14) 

(2.15) 


Since we have assumed that T\ and are two independent slow time scales in A{Ti, T 2 ), 
we have the following relation when the original time scale t is used as the temporal 
parameter for A{f) and B{f): 


r, A- 

Uf,-^ — V . .. 

'dTi 




dA 
, dB 


^ dB 

d,B = p—+p~ — 


(2.16) 

(2.17) 


We now combine the solvability conditions at ©(ry) and at 0{ri^) by using the above two 
relations, and write both equations in the original time units (also recall that / = ry/o 
and 7 = 770 ). For a frequency ratio m/n= 1/2, we have 


dfA = 


2ujq 


P f 3r^ 2(1 — r)^ 1 2 

— — +---r(l-r)cos(/) -7 


w. 


B 


rf r(l — r)p . , 

- 3 -sm (/) — 7 ) A, 


LOO 


2 a;, 


0 


dfB = -- 


2ojo 


P f 3r^ 2(1 — rY . 2 

— ( — H---h r(l - r) cos ^ - p 


Uln 


A 


"/ r(l-r)p 


OJo 


2iOQ 


sin(/ + 7 B. 


(2.18) 


(2.19) 


The eigenvalues of the above linear system determine the linear stability of the flat 
surface for subharmonic resonance. The condition for linear instability is that at least 
one eigenvalue has a positive real part. This condition reads. 



( 2 . 20 ) 


where we have grouped the dimensionless damping parameter ( 7 /a;o) and forcing ampli¬ 
tude {rf /up) together, which are actually the expansion parameters. The RHS of the 
above equation is of higher order in ^/ujq or in rf /loq than the LHS. However, as we 
shall see later (1 — r)/r « 4 for certain values of </> along the bicritical line; thus the 
first term on the RHS has a quite large coefficient: 37.0(r//a;§)‘*. Because of that, we 
do not neglect higher order terms on the RHS unless we consider the case of extremely 
weak damping j/ujo ^ 0.16. We have also checked that corrections to the threshold 
value of the driving force from orders higher than 0{p) do not contain terms of the 
form of ((1 — r)/r)"(r//a;§)" for n > 4. Therefore the truncation at 0{p) is a good 
approximation. 
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For frequency ratios min ^ 1/2, the instability condition reads, 




mujQ 


> 


2 f2 


r f 


2mr 


1 — r 


mojo 


■ ( 2 . 21 ) 


2 . 2 . 2 . Sq = 

In this case, the condition for the Faraday instability up to 0{r]) reads, 

(!-»')/ 


nojo 


> 7- 


( 2 . 22 ) 


The instability condition up to 0{rf') for the case of frequency ratio mjn = 1/2 \s given 

by, 


( 1 -?-)/ 

Au;l 


2r 
\ — r 


(I-’’)/ 




sinffl — 


> 


2 r 


3 2 

2 “ 3 V1-r 


7 

2 wo 
2 ' 


( 1 - 0 / 


(1 -^)/ y _ ( 


l — r 


4u;2 


V 2 wc 


(2.23) 


In this case, since r/(I — r) is expected to be small (in Muller’s experiment, r/(l — r) < 
0.45) along the bicritical line, the approximate threshold value of the driving force fth 
can be found with relatively good accuracy by solving Eq. (^.23) perturbatively, 


(1 - r)fth ^ _ 1 _ 

dWg 2ui0 


bi 


(— 

\ 2 wo 


(—)' 

\ 2 wo/ 


where 


(2.24) 


62 = 10 


l — r 


sin 9+2 


1 4 

4 3 V 1 — 


bi = -2 
2 


l — r 


-2 


sine/, 
4 


l—r 


In terms of the driving amplitude gz, the threshold value reads, 


I2h 16z/fcl2?iW0 
9z -I 

l — r 


l + h-^+b2 

LOo 


LOo 


where ki 2 h is determined by goki 2 h + +b^\ 2 hlP ~ ^^o- 

For frequency ratios that m/n 1/2, the instability condition reads, 


(i-o/ y / 7 y 

\nuio ) 


u'^LOq 


> 


{l-r)^f /3 


2n^ 




\1 — r 


(2.25) 

(2.26) 

(2.27) 


-(—)' 

\nu;o J 


The threshold value of the driving amplitude gz reads, 


(„) _ Smoknoop 


91 ' = 


l — r 


1 / 1 

1 + ^ U- 


2n^ 


8 \ 2 v? — mr \1 — r 


f 

V ncop J 


(2.28) 

(2.29) 


where kn is determined hy gpkn + Tk/^/p = n’^cop. 

For a frequency ratio m/n = 1/2, subharmonic resonance and harmonic resonance will 
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occur simultaneously (for the same value of the driving amplitude) when 
The value of is given by Eq. (2.20), which is written as follows, 

n 2 


^ ^ ' 1 2 ^ 
fs - —fsSm(j)--ff = - 



Is - It 


1 — r 
2 r 


fl (2.30) 


where fs = rgi^^®^fci 2 s/( 4 a;Q) and 7 s = 2 vk'f2sl‘^o- The value of is given by 

Eq. (|2.23), which can be written as. 


fl + 


2 r 

1 — r 


f^smcf-jh = ^ 


2 r '' 

3 V 1^, 


r 2 „,2 

/ h T/i 


-4 


rfh 

1 — r 


(2.31) 


where //i = (1 — r)gi^‘^^'^ki 2 h/ (16wq) and 7 ^ = The bicritical line rb((()) can be 

obtained directly by solving Eqs. ( |2.3C ) and (2.31) numerically. Alternatively, rf,((()) can 
be obtain ed by substituting the expr ession of of Eq. ( 2.27 ) as well as that of 
into Eq. ( 2.30[) , and then solve Eq. ( 2.30 ). Except for a small difference near cj) = 7 r/ 2 , 
where r is largest on the bicritical line, the second way is much easier. In the results 
presented below, the values of rb(0) directly calculated from Eq. ( 2.30| ) and ( ^.3l| ) are 
used. 

The solid curve in Figure is the bicritical line calculated by using parameters ap¬ 
propriate for the experimental data of Muller (1993)| . Note that for simplicity we have 
considered surface waves in the infinite depth limit in the above calculations. By assum¬ 
ing that the viscous dissipation near the bottom boundary layer can be neglected, our 
results can be generalized to the case of a fluid layer of finite depth d by replacing the in¬ 
finite depth dispersion relation by uj{k)'^ = taTih.{kd){gok + Tk^/p) and gz by gz tanh(fc(i). 
Since ki 2 sd = 2.0 and ki 2 hd = 3.9 in Muller’s experiments, such finite depth corrections 
are small, but have been included in Fig. ^ We also note that the value of ki 2 s calculated 
from the hnite depth dispersion relation using appropriate fluid and forcing parameters 
is different from that observed experimentally by Muller, and we have used the experi¬ 
mental value of wavenumber ki 2 s 0. Our analytical results agree qualitatively with the 
experimental results of Muller, shown as gray symbols. Other curves in Fig. ^ are bicrit¬ 
ical lines for the same values of Wq, fci 2 sj and ki 2 h as the solid curve, but v = 0.15 for 
the dotted curve, v = 0.1 for the dashed curve, and v = 0.05 for the long-dashed curve. 
Therefore our calculation is in qualitative agreement with the bicritical line determined 
experimentally. The small quantitative discrepancy between our analytical results and 
experimental results is probably due to the relatively large values of the damping pa¬ 
rameters ( 7 s = 0.17 and 7/1 = 0.33). It is likely that the quasi-potential approximation 
is not accurate in this parameter range. 

For frequency ratios m/n 1/2, we predict that the bicritical line is independent of 
the phase difference (j), and can be readily obtained from Eq. ( ^.21| ) and ( ^.28| ). Further 
experimental studies of bicritical lines for different values of frequency ratio, v, ujq, and 
surface tension T would be interesting to provide additional tests of our results. 


f The value of ki 2 s calculated from the finite depth dispersion relation using appropriate fluid 
and forcing parameters (p — 0.95 > g/cm®, F = 20.6dyn/cm, d — 0.23cm, and wo = 27rx27.9Hz) 
is lOcrn”^, or Ai 2 s = 0.63cm. This value of the critical wavelength is significantly different from 
what was observed by Muller, Ai 2 s ~ 0.72 cm. The calculated value of k\ 2 h does agree well with 
the observed value, ki 2 h = 17cm“^, or \\ 2 h ~ 0.37cm. 
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3. Standing Wave Amplitude Equations 

We present in this section a weakly nonlinear analysis of parametric surface waves 
driven by two frequencies, extending the calculations presented earlier for the single 
frequency case by Zhang fc Vinals (1996b)| . For fluids of low viscosity, the equations 
governing fluid flow and the boundary conditions at the free surface were expanded in 
the (small) width of the vortical layer adjacent to the free surface. The bulk flow is 
then potential, but satisfies effective boundary conditions on the moving surface. We 
furthermore neglected viscous terms that are nonlinear in the free surface variables to 
arrive at the so-called Linear Damping Quasi-Potential equations. They involve only the 
surface’s deviation away from planarity h, and the surface velocity potential <i>, but no 
longer depend on the bulk velocity field, 


dth{x, t) = -/y h + V<^?-V ■ (W$) -t- -y^{h^V<S>) 


-V{hV^)+V 


hT>{hV‘^>) + 


(3.1) 


dt$(x, t) = 7 V^< 1 > — (Go — FqV^) h — if ^sin 2 t -|- Ksm{2pt + 0 )^ h 


_ 1 (V$)2 - (f)$) hy^<i> + t>{h'D<i>) - y V • (V/i(V/i)2). (3.2) 


Linearization of these two equations leads to Eq. ( ^.21) . When the forcing component 
sin 2 ma'of^ dominates, we choose l/(ma;o) as the unit of time and 1/kom as the unit of 
length with gokom + = mfojQ- For the other case when the forcing component 

sui2nLOot dominates, we choose l/(nwo) as the unit of time and 1/kon as the unit of 
length with gokon + = u^ujq. 2 ? is a linear and nonlocal operator that multiplies 

each Fourier component of a field by its wavenumber modulus (Zhang & Vinals 1996b). 
We can write the system of equations for the two cases in the same dimensionless form, 
with the dimensionless variables 7 , Go, Fq, /, k, and p having different values for the 
two cases. When the forcing component siii2muJot dominates, we have 


7 = 




Gn = 


mujQ 

2 2 ’ 


/ = 


9z ^Om'^ 


1 — r n 

K = -, p = —, 


ro = 


1 /in. 


Om 


9 2 ’ 
pm^ujQ 


(Go -f Fq — 1), 0 — 4 >. 


For the other case (the forcing component sin2na;ot dominates), we have 


7 = 

Go = 


nujo 

gokon 


f = 


gzkon{l - r) 




-j-, P= — 

1 — r n 


Fn = 


Tkl 


On 


9 2 ’ ■*■0 9 2 ’ 

pn^ujQ 


(Go -|- Fo — 1), 


0 = - 


m 


(3.3) 

(3.4) 

(3.5) 

(3.6) 


3.1. Derivation of SWAE’s 

As mentioned above, the detailed procedure for the derivation of standing wave amplitude 
equations parallels that presented by Zhang fc Vinals (1996b)| for the case of sinusoidal 
forcing. We seek nonlinear standing wave solutions of Faraday waves near onset, i.e. 
e = [f — j)h ^ 1. The quasi-potential equations (Eqs. Q and (|3.2D) are expanded 
consistently in with multiple time scales, 

/i(x, t, T) = e^/^/ii(x, t, T) + e/i 2 + -|-, 


(3.7) 
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$(x, t, T) = (x, t, T) + e $2 + + • ■ •, 
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where hi and $i are the linear neutral solutions and can be found in a similar way to 
the case of sinusoidal forcing. For simplicity, we shall only consider terms up to order / 
or 7 for hi and $ 1 . They are 


hi(x, t) = cost + ^ sin3t + ^sin^( 2 p + l)t + 0 j 

+ ^y^^^ 3 ^sin^( 2 p-l)t + 0j ^ (T) exp • x) + c.c , (3.9) 

$i(x, t) = — sint + 7 cost + — cos3t + ^Qg j + l)t + 0 ) 

4 2p{p+l) \ ) 


7 k( 2 p - 1 ) 

2p{p - 1) 


cos^(2p — l)t + 0^1 ^ I^Aj(r) exp • x^ + c.c , (3.10) 

' ' ' 1=1 


where T = et. 

At 0(e), we have the following non-homogeneous linear equation for / 12 , 


dtth 2 — 27V^9t^2 + (Go — Tq V^)T>/i2 + 47 ^sin 2t + k sm{2pt + 0)^ 2?/i2 

N e _ _ _ 

= ^ y^2(l + Cji) — cos2t 1 + Cji y^2(l + Cji) 

—7sin2t (5/2 + cji) ^1 + cji — ^2(1 + cji)^ H- -^-^^2(1 + cji) 

+ Y~^ sin(2pt + 0) —^2(1 + Cji) + ^1 + Cji - ^2(1 + Cji)^ 

~~ sin[2(p+l)t + 0] ^ {l+Cji)^2{l+Cji) + (p+1) (^ + Cji — ^2{\+Cji^ 

+ ~ sin[2(p—l)t+0] ———(1 + Cj7)-^2(1+Cj/) + (p—1) ^l+c^; — •^2(1 + Cj7)^ 

X AjAi exp {iikj + ki) ■ x^ + c.c. 

v r _ P _ 

+ 1-cp --^^\/2(1^-^ 

I 

-7sin2t (5/2 - cp) ^1 - cp - 1^2(1-cp)^ + ^ 


+ T~^ sin(2pt + 0) 
1 — 


2(1 Cji) -\- p (1 Cj7 \/^(^ 


-^sin[2(p+l)t + 0] ^_||~^^ (l-Cj-;)y2(l-cp) + (p+1) ^2{l-Cji^ 

+ ^sin[2(p-l)t + 0] ‘^:^^{l-Cji)^2{l-Cji) + (p-1) ^l-cp-y2(l-cp)^ | 
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X AjA* exp (^i{kj — ki) ■ 


+ c.c. 


.(3.11) 


Before obtaining the detailed form of the equations, let us mention that the generic 
form of the SWAE’s for two-frequency driven Faraday waves can be obtained by symme¬ 
try considerations ( Edwards fc Fauve 1994). For the case of single frequency forcing the 
SWAE’s derived in Zhang fc Vinals (1996b) contain only third order nonlinear terms, 
but no quadratic terms. The exclusion of quadratic terms can be understood there from 
the requirement of sign invariance of the SWAEs. Subharmonic response of the fluid 
surface to the driving force /sin(2u;ot) implies /ij(x,t -|- tt/loq) = Here hj is a 

linear unstable mode, 

hj = ^coswot -I- ^ sin Swot + • ■ 0^ + c-c.^ , (3.12) 

where only odd multiples of frequency wq appear. As a result, a sign change of Aj 
is equivalent to a time displacement in a period of the driving force, t —>■ t + tt/loq. 
Because of the invariance of the surface wave system under such a time displacement, 
the amplitude equation of Aj must be sign invariant, which obviously excludes quadratic 
terms. However, quadratic terms can arise with two-frequency forcing. If the frequency 
ratio m/n = even/odd or odd/even, the SWAE’s are sign invariant if the odd frequency 
dominates, and otherwise if the even frequency dominates. When the frequency ratio 
m/n = odd/odd, the SWAE’s are sign invariant. A general consequence of the loss of 
sign invariance is the appearance of quadratic terms in the amplitude equation. 

We concentrate here on the case of frequency ratio m/n = 1/2, since detailed exper¬ 
imental results are available for this case. Similar calculations can be done for other 
frequency ratios, and they might be interesting for future experimental studies. We 
further restrict our attention to the case that the odd frequency dominates, which also 
corresponds to the subharmonic side of the bicritical line. In this case, the solutions at 
0(e^/^) are. 


hi = 


cost -|- ^ sin3t -I- ^ sin(3t -I- /))j ^ \Aj exp {kj -x^ -l-c.c , (3.13) 


$1 = 


37 37K 


1=1 


— sin t -I- 7 cos t -I —7- cos3t H -— coslSt + d)) 

4 4 ' 


\Aj exp {kj -x^ -l-c.c ,(3.14) 


1=1 


where k = (1 — r)/r. At 0{e), we have, 


dtth 2 - 27V^9t/i2 + (Go - roV^)T>/i2 + 47^sin2t -|- «;sin(4t -f (/) ]Vh 2 


jsf r 

= ^ ^ +Cji)-cos2t 


3,1=-^ 


—7 sin 2t 


1 -I- Cjl - 


3 Cjl 


\J2{1 Cjl) 


—— sin(/ ^1 -\- Cjl (2 -\- 'icji)-\j2{\ -\- Cji^ 
(5/2 -I- Cjl) ^1 -I- Cjl - ^2(1 -I- Cji)'^ -\ ^^\/2(l + Cjl) 


— — cos (j) [ 1 + Cjl + {2 + 3cji) 


2(1 -I- Cjl) 


AjAi exp (^i{kj -I- fc;) • x^ 


-I- c.c. 

























n 

--cos 
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i,J=i 


sin ^ 1^1 - Cj/ + (2 - 3cji)^2(l - Cji)^ 
(5/2 - c,0 (^1 - c,i - ^2(1 - c,0) + ^^72(1 -c,0 


—7 sin 21 


(1 _ + (2 - 3c,)/iirT7;i) ] I »p (.(i, - i,) ^ x) + „]3.13) 


We now solve the above equation for /i 2 - Since there are no secular terms on the 
RHS, we are only interested in the particular solution for /12 caused by the RHS of that 
equation. The particular solution reads, 


JV X 

^2 = ^ -I (aji + fdji cos 21 + jSji sin 21 ) AjAi exp (^i{kj + ki) ■ + c.c. 

j,i=i ^ 

+ + Pji cos 21 + ^5ji sin 2l) AjA* exp (^{kj — fc;) • + c.c. |, (3.16) 


where 


1 + Cjl 


27^ 




-H 


Pji = 


^ 4 [Go + 2ro(l + Cj/)] Go + 2ro(l + Cjj) 

'ji {^oi - + 7 ^ (8{1 + Cjl) - 2 ka / 2(1 + cji) cos(j)j Nji 

72 ^8(1 + Cjl) - 2ka/2(1 + Cji) cos + E'ji - 8 ^‘^EjiMji 


Sji = -- 


^ 8(1 + Cjl) — 2 s; 1 ^ 2(1 + Cjl) cos<f>^ Hji + EjiNji 
72 ^ 8(1 + Cji) - 2 ka/ 2(1 + Cjl) cos (j)j + E^i - Sj'^EjiMji 


Mu = 


\/2(l + Cjl) 


(3.17) 

(3.18) 

(3.19) 

(3.20) 


^ Gq H- 2ro(l + Cjl) 

Dji = [Go + 2ro(l + Cjl)] ^2(1 + Cji) - 4, (3.21) 


Eji = Dji + 2 kj^J 2{1 + Cjl) sin(/), 


(3.22) 


Hp = 1 
and 


+ Cj; y^2(1 + Cjl) ^ + (2 + 3 cj 7 )y^ 2 (l + Cji)^ sincj), (3.23) 


^ji “ ^2 ~ ^2(1 + Cji)^ H -^-^^2{1 + Cjl) + (1 + Cji)Mji 

— — ^1 + Cjl + (2 + 'icji)^J 2(1 + Cji^ cos(/.(3.24) 


aji, (dji, and 6ji can be obtained by replacing Cji with —Cji in the expressions for aji, 
(dji, and Sji respectively. 
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The solution for 4)2 is, 


iV ^ 

4)2 = •! {jUji + jVji cos2t + Wji sin2t) AjAi exp (^i{kj + ki) ■ + c.c. 

+ cos2t + Wji sin2t) AjA* exp (^i{kj — ki) ■ + c.c. (3.25) 


where 


— 2 ■*" ~ Y^2(1 + Cji), 

”11 = J + (ft. -1) ^2(1+ '..)+ ^2(l’ic„y 


1 1 


2/3,7 




(3.26) 

(3.27) 

(3.28) 


and Uji, Vji, and Wji can be obtained by replacing Cji with —Cji in the expressions for 
Uji, Vji, and Wji respectively. 

We note that the triad resonant condition is Eji = 0 in Eq. ( ^.22| ), which is different 
from that of the case of sinusoidal forcing {Dji = 0) (Zhang & Vinals 1996b). This 
difference is important since it results in different values of Cji at the triad resonance, 
which is now also a function of the damping parameter 7 , the ratio of amplitudes r, 
and the phase difference (p. As an example. Figure || shows the modified triad resonant 
condition for Tq = 1, 7 = 0.15, and r = 0.3. As we see that the value of is 
close to 90° for some values of (p. From the results in the case of sinusoidal forcing, 
we know that when the wavevectors of two standing waves are separated by an angle 

(r) 

of 0 };^ the pattern is strongly suppressed. Two waves of wavevectors separated by this 
angle will excite a linearly stable mode with an amplitude inversely proportional to the 
damping coefficient. Energy is then dissipated by the stable mode, and hence the two 
original waves are effectively damped as compared to other unstable modes that do not 
satisfy the triad resonance condition. Thus, without further calculations one already can 
conclude that square patterns are not favored when the phase difference p is close to 7r/2 
(or if close to 0 in Miiller’s notation) for this set of parameters. 

At the calculation is again similar to the case of single frequency forcing. The 

solvability condition gives the following standing wave amplitude equations 


\^-A- 
j dT ~ ^ 


5(1)1^,f 


N 




where j = 1 , 2 ,' 


,/V, 


, , 28 + 9ro + (12 + 9ro)KCOS(() „ 3(1 + kcos(()) 

5(1) =-271-+ 2a,-, +--- 


64 


/3,7 o ! 


(3.29) 


(3.30) 


and 


^ 3r,(a.cos0) ( 3 - 


32 


+ ^1 + Cji - ^2(1 + Cji)^ ^ 


+ 1 Cji V 2(1 Cji) 


1 + K cos P 
4 

1 + K COS 


-Wji - Vji 

i _ 

-Wji - Vjl 
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/I 3(1 + KCOS(/l) ^ 1- \ 

+(1 + Cji) l^aji H-g- Pji - 2^jiJ 

N 3(1 + Kcosth) - 1 7 = \ N 

+ ~ '^ji) “I ^ t^ji ~ 

When K = 0, Eq. ( |1.29| ) reduces to the standing wave amplitude equation for sinusoidal 
forcing as expected (Zhang & Vihals 1996b). 

In contrast to single frequency forced Faraday waves, we note that 5 ( 1 ) < 0 for some 
ranges of parameter values in this two frequency case. When 5 ( 1 ) < 0, the amplitude 
Aj would increase without limit. Therefore higher order terms (at least of fifth order) 
will be required in the SWAE to saturate the amplitude. The regions where g{l) < 0 for 
particular sets of parameters are shown as black in Fig. ^ The steady state of parametric 
surface waves in parameter regions such that (/(I) < 0 cannot be determined from the 
SWAEs (Eq. ( |3.29D ). Analytical calculations of the relevant higher order terms in this 
case are algebraically much more complicated than the calculations presented here, and 
we have not done such calculations. We shall restrict ourselves to the parameter range 

that gives a positive ( 7 ( 1 ) in further analysis of this section. _ 

When g{l) > 0, we can rescale the amplitude Aj in Eq. ( ^3.291 ) as Aj = g{l)Aj. We 
have the following standing wave amplitude equation for the scaled amplitude. 


7 dT ^ 


|AjP+ ^ 5(cj/)|A/p 


A, 


(3.32) 


where g{cji) = g{cji)/g{l). Note that g{cji ±1) = 2. 

The scaled nonlinear interaction coefficient g(cji) (we have suppressed the tilde since 
we will only refer to the scaled nonlinear coefficient in what follows) for the regions with 
a positive 5 ( 1 ) depends strongly on the phase difference (j). This is due to the effect of 
the modified triad resonant condition (see Fig. ||), as well as other effects, such as the 
dependence of the surface wave amplitude at the modified triad resonance on (j). As 
an example, Fig. ||(a) shows the scaled coefficient g{cji) as a function of Cji for purely 
capillary waves with 7 = 0.1 and r = 0.25. The general trend in these curves can 
be understood from the different triad resonant condition for different values of (j). As 
4> decreases from 37r/2 to 7r/2, the angle (recall that Cji = cos(9ji)) at the triad 
resonance increases to a value close to 7r/2 (see Fig. ||), and therefore the values of 
g{cji) changes from a minimum (< 1) to a local maximum (> 1). This change in gicji) 
makes a pattern of square symmetry unstable for (f close to 7r/2 (since g(0) > 1). By 
increasing the values of r in Fig. ib), the amplitude of the smaller forcing component 
(Ksin(4t + (f)) decreases, and therefore the changes in g{cji) for different values of </> are 
smaller. We note that with the damping parameters used in Fig. |^, surface waves may 
have synchronous (harmonic) response to the driving force, i.e. below the bicritical line 
in Fig. 1^ for some values of the phase difference </>. In that case, the above results of 
g{cji) become irrelevant for these values of ip. 

For larger values of the damping parameter 7 , triad resonant interactions are more 
strongly damped, as is the case for single frequency forcing. However, the situation is 
more complicated here since the modified triad resonant condition depends on 7 (actually 
it depends on the driving amplitude /, which is approximately equal to 7 at onset). For 
example, larger values of 7 can make the angle at the triad resonance closer to 7r/2. 
Figure shows the function g{cji) for the same values of parameters as in Fig. ^ except 
the value of 7 is raised to 7 = 0.15. 
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3.2. Pattern selection near onset 


A qualitative analysis of pattern selection for the two frequency case was already advanced 
by Muller (1993), based on general symmetry considerations and a typical shape of the 
nonlinear coupling function g{cji). The derivation of this function given in the previous 
section allows us to obtain quantitative predictions for the regions in parameter space 
in which patterns of a given symmetry minimize the appropriate Lyapunov function for 
this problem. Since Eq. ( |3.32 ) is of gradient form Ij^dxAj = —dtFIdA*, a Lyapunov 
function J- can be defined as, 


N 


i=i 




N 




Since 


dE 

dr 


N 

= E 

i=i 


dE 

dA, 


dTA-i 


dE 

dA* 


N 


dr A* = < 0, 


(3.33) 


(3.34) 


1=1 


the only possible limiting cases of such a dissipative system, in the limit T —> oo, are 
stationary states for the amplitudes Aj. Only the states which correspond to local minima 

of the Lyapunov function are linearly stable. _ 

Apart from the trivial solution of Aj = 0 for j = 1, • • •, A^, Eq. ( 3.32 ) has a family of 
stationary solutions differing in the total number of standing waves N for which Aj ^ 0. 
By considering t he ca se in which the magnitudes of all standing waves are the same, i.e. 
|Aj| = jAj, Eq. ( 3.32| ) has the following solutions, 


N 


- 1/2 


l^jl — 1^1 — I ^ E 9{cji) 

V 1=RIJ3) 

The values of the Lyapunov function for these solutions are. 




9{cji) 


(3.35) 


(3.36) 


We shall only consider pattern structures for which the angle between any two adja¬ 
cent wavevectors kj and kj+i is the same and amounts to tt/N (regular pattern). We 
summarize our results concerning regular patterns in Fig. Note that our results only 
apply to the subharmonic response to the driving force, which corresponds to the side of 
bicritical line with smaller amplitude ratio r. Different regions are labeled by the pattern 
structure that has the lowest value of the Lyapunov function En except the region shown 
in black, in which the self-interaction coefficient (/(I) is negative. Fig. Sa) uses all the 
appropriate fluid and forcing parameters from Muller’s experiment (Muller 1993). As 
shown in Fig. |(a), hexagonal or triangular patterns have the lowest value of Lyapunov 
function in a region around (j) = nj2 and close to the bicritical line while square pat¬ 
tern is stable near the bicritical line and around cj) = 37r/2 (or —t: 12). This result is in 
qualitative agreement with Muller’s experiment although our theory predicts a smaller 
stable hexagonal/triangular region for this set of parameter values. In a region close to 
the bicritical line and with 0 < 0 < 7r/2, differences in the values of Lyapunov function 
for different patterns become smaller, and most stable patterns of different symmetries 
(including quasicrystalline ones) correspond to neighboring smaller parameter regions. 
It is possible that the preferred pattern may not be selected due to such small differences 
in Lyapunov function for patterns of different symmetries, or there is no selection at all. 
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This result is in partial agreement with Muller’s experiment, where he found disordered 
patterns in a parameter region overlapping with the above-mentioned region. 

Fig. ||(b-d) uses the same parameter values except the fluid viscosity which is changed 
from V = 0.20cm^/s to 0.15cm^/s, O.lOcm^/s, and 0.05cm^/s respectively. We note that 
as V decreases (smaller viscous damping), the hexagonal/triangular region becomes larger 
and the center of this region is also shifted towards a larger value of (f. For v = 0.05cm^/s, 
the hexagonal/triangular region reaches larger values of r > 0.60. Currently there are no 
experimental results known to us that can be used to test our predictions as the damping 
parameter varies. 


4. Discussion and snmmary 

We comment further on the approximations considered to obtain our main results 
presented in Fig. In the quasi-potential approximation, terms of the order of 7 ^/^ or 
higher have been neglected. Consistent with that, we use approximate linear solutions, 
correct up to order / or 7 (the perturbation expansion for the linear solutions) in deriving 
the SWAEs (Eq. (3.29)). On the other hand, we included terms of the order of p, /q. 


or p in the linear stability analysis to obtain the bicritical line. We considered these 
terms because of the large coefficients for such terms in the perturbation expansion for 
the bicritical line in order to compare our results with Muller’s experiment which used 
relatively large values of 7 (especially for the case of harmonic response = 0.33). A 
fully consistent calculation correct up to the order of p will be much more difficult. 
Based on the qualitative agreement of the bicritical line we obtained with experiments, 
we expect that a fully consistent calculation would change our results quantitatively, but 
not qualitatively. 

In summary, analytical results of bicritical lines rb{(f>) are obtained by using a multiscale 
perturbation expansion. The results of rtp) for frequency ratio 1/2 agree qualitatively 
with Muller’s experimental results. We also derive the SWAEs for frequency ratio 1/2 
(for the case of subharmonic response). We found that the triad resonance condition 
is modified due the presence of the second frequency component. For certain values of 

(r) 

the relative phase difference between the two forcing components, we found that 0 ); 
becomes close to 90°. As a result, square patterns become unstable in this parameter 
region. Even though quadratic terms are prohibited for subharmonic responses, hexag¬ 
onal or triangular patterns can be stabilized with the presence of the second frequency 
component, which is in agreement with experiments. We also studied pattern selection 
for different values of the damping parameters, and found that hexagonal/triangular pat¬ 
terns are stabilized in a larger region for smaller values of the damping parameter. This 
is a prediction that awaits experimental verification. 

This research has been supported by the U.S. Department of Energy, contract No. 
DE-FG05-95ER14566, and also in part by the Supercomputer Computations Research 
Institute, which is partially funded by the U.S. Department of Energy, contract No. 
DE-FC05-85ER25000. 
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Figure 1. Bicritical lines of Faraday waves driven by a two-frequency force of frequency ratio 
m/n = 1/2. The gray symbols are taken from Muller’s experiment (Muller 1993), and the solid 
curve is calculated using parameters of Muller’s experiment: v = 0.2cm^/s, wqI’I'k = 27.9Hz, 
ki 2 a = 8.7cm“'^, and ki 2 h = 17.0cm“^. Other curves are calculated using the same values of loq, 
ki 2 a, and ki 2 h as the solid curve, but v = 0.15cm^/s for the dotted curve, v = O.lcm^/s for the 
dashed curve, and v = 0.05cm^/s for the long-dashed curve. 
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TT 



4 ^ 


Figure 2. Modified triad resonant condition, as a function of the phase difference cj), for 
the subharmonic response in two-frequency Faraday waves of frequency ratio 1:2 with Fq = 1, 
7 = 0.15, and r = 0.3. The dashed line corresponds to the value of for purely capillary 
waves in the case of sinusoidal forcing. 
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Figure 3. Nonlinear coefficient gicji) of the SWAEs for two-frequency driven Faraday waves 
for purely capillary waves with damping parameter 7 = 0.1 and the relative amplitude r = 0.25 
in (a), r = 0.35 in (b). The phase difference f = 'k 12 for the thick solid curves, f = O.Gtt for the 
thick dotted curves, f = O.Tn for the thick dashed curves, f = O.Stt for the thick long-dashed 
curves, (j> = O.Ott for the thick dot-dashed curves, (j> = -k for the thin solid curves, (j> = 1 . 37 r for 
the thin dotted curves, and (^ = 37 r /2 for the thin dashed curves. 
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Figure 4. Nonlinear coefficient g{cji) of the SWAEs for two-frequency driven Faraday waves 
for purely capillary waves with damping parameter 7 = 0.15 and the relative amplitude r = 0.25 
in (a), r = 0.35 in (b). The phase difference (j> = ^/2 for the thick solid curves, 4> = O.Ott for the 
thick dotted curves, (j> = O.Ttt for the thick dashed curves, (j> = O.Stt for the thick long-dashed 
curves, (j> = O.Ott for the thick dot-dashed curves, 0 = tt for the thin solid curves, (j> = I.Stt for 
the thin dotted curves, and = J,'k for the thin dashed curves. 
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Figure 5. Regions labeled by the pattern structure that has the lowest value of the Lyapunov 
function: L - lines, S - square, H/T - hexagon/triangle, and Q - quasicrystalline patterns 
(8-fold, 10-fold, 12-fold, ...). The black region indicates where the self-interaction coefficient g{l) 
is negative. Figure (a) is calculated using parameters of Muller’s experiment: p = 0.95 g/cm®, 
F = 20.6 dyn/cm, v = 0.20 cm^/s, wo/Stt = 27.9 Hz, 2-k jkvia = 0.72cm, and 2'iT/ki2h = 17.0 cm. 
Figure (b), (c), and (d) are calculated using the same values of p, F, ojq, fei 2 s, and as (a) 
but V = 0.15cm^/s, v = O.lcm^/s, and v = 0.05cm^/s respectively. Note the different r scale 
used for (d). 












